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ABSTRACT 

The neutron component is likely an inevitable ingredient of a gamma-ray burst (GRB) baryonic 
fireball, in essentially all progenitor scenarios. The suggestion that the neutron composition may alter 
the early afterglow behavior has been proposed, but there is no detailed calculation so far. In this 
paper, within the popular internal shock scenario of GRBs, we calculate the early optical afterglow 
lightcurves of a neutron-fed GRB fireball for different assumed neutron fractions in the fireball and 
for both ISM- and wind-interaction models. The cases for both long and short GRBs are considered. 
We show that as long as the neutron fraction is significant (e.g. the number of neutrons is comparable 
to that of protons), rich afterglow signatures would show up. For a constant density (ISM) model, 
a neutron-rich early afterglow is characterized by a slowly rising lightcurve followed by a sharp re- 
brightening bump caused by collision between the leading neutron decay trail ejecta and the trailing 
ion ejecta. For a massive star stellar-wind model, the neutron-rich early afterglow shows an extended 
plateau lasting for about 100 seconds before the lightcurve starts to decay. The plateau is mainly 
attributed to the emission from the unshocked neutron decay trail. When the overlapping of the initial 
prompt 7— rays with the shocks and the trail is important, as is common for the wind model and is 
also possible in the ISM model under some conditions, the IC cooling effect suppresses the very early 
optical afterglow significantly, making the neutron-fed signature dimmer. For short GRBs powered 
by compact star mergers, a neutron-decay-induced step-like re-brightening is predicted, although the 
amplitude is not large. All these neutron-fed signatures are likely detectable by the Ultraviolet Optical 
Telescope (UVOT) on board the Swift observatory if GRB fireballs are indeed baryonic and neutron- 
rich. Close monitoring of early afterglows from 10s to 1000s of seconds, when combined with detailed 
theoretical modeling, could be used to potentially diagnose the existence of the neutron component 
in GRB fireballs. 

Subject headings: gamma rays: bursts — radiation mechanisms: non-thermal — shock waves 



1. INTRODUCTION 

The suggestion that Gamma-ray burst (GRB) fire- 
balls should contain a good fraction of neutrons has at- 
tracted broad attention recently, since in essentially all 
progenitor scenarios the neutron component is likely an 
inevitable ingredient for a baryonic GRB fireball (e.g., 
Derishev et al. 1999a; Beloborodov 2003a; Pruet et 
al. 2003). For instance, core-collapse of massive stars 
would lead to an outflow from an iron-rich core with the 
parameter x, the neutron-to-proton number ratio, be- 
ing > 1 (e.g., Beloborodov 2003a). In the neutron star 
merger model that may be valid for short-hard GRBs, 
one would also expect x > 1. Photo dissociation during 
collapse or merger, as well as n, p decoupling and inelas- 
tic collisions, would both drive x towards unity, although 
such an equalization process is likely to remain incom- 
plete (Bachall & Meszaros 2000). Weak interactions in- 
duced by the intense neutrino flux from the central engine 
can result in significant proton-to-neutron conversion, es- 
pecially if resonant neutrino flavor transformation takes 
place (Qian et al. 1993; Fuller et al. 2000). 

Derishev et al. (1999a, 1999b) first investigated the dy- 
namics and the possible observational signatures of a rel- 
ativistic neutron-rich fireball. This was followed by many 
related investigations. One advantage of the neutron rich 
model is that the baryon-loading problem for GRBs can 
be ameliorated if significant fraction of baryons confined 



in the fireball are converted to neutrons (Fuller et al. 
2000). The existence of the neutron component likely 
leaves various observational signatures. For example, the 
decoupling between the neutron and the proton compo- 
nents during the early fireball acceleration phase would 
give rise to a distinct multi-GeV neutrino emission sig- 
nature because of the inelastic neutron-ion collisions in 
the fireball (Bachall & Meszaros 2000; Meszaros & Rees 
2000). Such a GeV neutrino signature is however not 
easy to detect in the near future. Lately Fan & Wei 
(2004) suggested that there should be a bright ultraviolet 
(UV) flash accompanying a GRB from neutron-rich in- 
ternal shocks, as long as the burst is long enough^. Such 
a UV flash is in principle detectable by the Ultraviolet- 
Optical Telescope (UVOT) on board the Swift observa- 
tory. However, since such signals happen early (typically 
> 14(1 -|- z)s after the burst trigger), given the nominal 

^ As reahzed in Fan et al. (2005a), the accelerated electrons are 
mainly cooled by the inverse Compton scattering with the initial 
prompt 7— ray photons (Beloborodov 2005), so that the UV flash 
may be dimmer by 5 mag, but it is still bright enough to be de- 
tected by the Swift UVOT. The predicted such low energy photon 
flash accompanying the prompt 7— ray emission may have already 
been detected in GRB 041219a (Blake et al. 2005; Vestrand et al. 
2005), for which it is found the optical emission during the burst 
is variable and correlated with the prompt 7— rays, indicating a 
common origin for the optical light and the 7— rays (Vestrand et 
al. 2005). This viewpoint has also been supported by the IR band 
observation (Blake et al. 2005). 
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UVOT on-target time (60-100 s), testing this signature 
requires an optimized detector configuration. The most 
promising approach to test the neutron component would 
be the signatures in the early afterglow phase, typically 
tens to hundreds of seconds after the burst trigger, when 
UVOT is likely on-target. Such early afterglow neutron 
signatures have been suggested earlier (Derishev et al. 
1999b; Beloborodov 2003b, hereafter BOS). However, 
there has been no detailed calculations so far, and this is 
the main focus of the current paper. 

Our work differs from BOS in three main aspects: (i) 
In BOS, the neutron shell (N-ejecta) always keeps ahead 
of the ion-ejecta (I-ejecta)^, so almost all the decayed 
products from the N-ejecta are deposited onto the ex- 
ternal medium and are used to accelerate the medium. 
Such an approximation is only valid when the burst du- 
ration is short enough (e.g. for short bursts). In this 
work we show that for typical long GRBs, the separation 
between the N-ejecta and the Tejecta is not clean. The 
N-ejecta would partially overlap the I-ejecta until a long 
distance R ~ several x R/3, where Rp = 8 x 10^^rn_2.5cm 
is the mean /?— decay {n ~* p + + Dc) radius of the 
N-ejecta, is the bulk Lorentz factor (LF) of the N- 
ejecta, and the convention — is adopted in cgs 

units here and throughout the text. Consequently, a sig- 
nificant fraction of the /3— decay products are deposited 
in the I-ejecta rather than in the external medium. So 
the neutron-decay trail (i.e., the external medium mixed 
with the /3— decay products) is less energetic than that 
suggested in BOS. Nonetheless, we confirmed the con- 
clusion of BOS that the presence of the neutron ejecta 
qualitatively changes the afterglow emission properties, 
(ii) BOS mainly discussed the dynamical evolution of a 
neutron-fed fireball. The energy dissipation rates in the 
neutron front and in the shock front have been calcu- 
lated. Although they can delineate the main character- 
istics of a neutron-fed fireball, these results are not easy 
to be directly compared with the observations. In this 
work, we extend BOS's investigations to calculate syn- 
chrotron radiation both from the shocks and from the 
neutron decay trail. We calculated in detail some sam- 
ple early optical lightcurves, which can be directly com- 
pared with the future observations by Swift UVOT or 
similar telescopes, (iii) BOS mainly discussed the influ- 
ence of neutrons on the forward shock (FS) emission. In 
this work besides the FS component, we also explicitly 
discuss the emission from the reverse shock (RS) region, 
which has been widely accepted to play an important 
role in shaping early afterglow lightcurves. 

The structure of the paper is as follows: In ^ we 
present the physical picture, including the neutron-rich 
internal shocks, the formation of the I-ejecta, the neutron 
decay trail and the dynamics of the neutron-fed fireball. 
In 50 (for long GRBs) and S|(for short GRBs), we model 
the dynamics of neutron-rich systems numerically and 
calculate the synchrotron radiation of various emission 
components. Sample early optical afterglows for typical 
parameters are calculated and presented. For long GRBs 
the cases for both a constant-density medium (ISM) 
and for a stellar-wind medium (wind) are presented. Our 

^ In most cases, ions are dissociated into protons. Following the 
convention in the literature, we call the proton ejecta also as ion 
ejecta. 



results are summarized in |5lwith some discussions. 

2. THE PHYSICAL PICTURE 
2.1. n-p decoupling in a neutron-rich fireball 

In the standard GRB model, a fireball made up of 7,e^ 
and an admixture of baryons is generated by the release 
of a large amount of energy i^iso > lO^^erg in a region of 
ro ~ lO^cm, where Eiso is the total energy of the burst 
assuming isotropic energy distribution. Data from the 
bursts with known redshifts indicate that a typical fire- 
ball is characterized by a wind luminosity L ~ 5 x lO^'^erg 
s^^ and a duration Tgo ^ 50s, measured in the observer 
frame. Above the fireball injection radius tq the bulk 
Lorentz factor (hereafter LF) F varies as F ~ r/rp ini- 
tially and saturates when it reaches an asymptotic value 
Tf < rj ^ constant, where rj is defined as 77 = L/Mc^ 
(e.g. Meszaros et al. 199S; Piran et al. 199S). 

For an n, p fireball, the picture becomes more com- 
plicated"^. The n and p components are cold in the 
co-moving frame, and remain well coupled until the co- 
moving nuclear elastic scattering time tjjp ~ (UpCroc)"^ 
becomes longer than the co-moving expansion time 
^cxp '"/cr, where ^ 3 x 10~^^ cm^ is the pion pro- 
duction cross section above the threshold ^ 140 MeV, 
rip is the co-moving proton density which according to 
mass conservation reads n'^ = L/[{1 + x)47rr^mpC'^r77]. 
Therefore the n, p decoupling occurs in the coasting 
or accelerating regimes depending on whether the di- 
mensionless entropy rj is below or above the critical 
value 77^ ~ S.9 x 10^Ll(^\~y'^[{l + x)/2]-i/4 (Bachall 
& Meszaros 2000; Beloborodov 200Sa). 

For r] < rjT^, both n and p coast with F ~ 77 = constant. 
For r] > rjj^, the condition t'^^ > t'^^^ is achieved at a 

radius 7'np/ro = rjTj{rij, /rj)^^/'^ . While protons are still 
being accelerated as Fp ~ r/ro, the neutrons are no 
longer accelerated and they coast at a LF of F ~ Ff_n ~ 

2.2 X 10^Lli^r~^^'^[{l + x)/2]~^^^V3^^^ (Derishev et al. 
1999a). From energy conservation, one gets the asymp- 
totic proton Lorentz factor Ff.p ~ 77(1 -|- x)[l — (x/[l + 
x])(6/7)(77^/77)^/^] (BachaU & Meszaros 2000). 

2.2. Internal shocks and formation of the I-ejecta 

In the standard fireball model, the long, complex GRBs 
are powered by the interaction of proton shells with dif- 
ferent LFs. The practical LF distribution of these shells 
could be very complicated, and detailed simulations of 
these interactions is beyond the scope of this work. As 
a simple toy model, here we follow Fan & Wei (2004) to 
assume that LFs of the shells follow an approximate bi- 
modal distribution (Guetta et al. 2001) with the typical 
values taken as rjf (fast shells) or t^s (slow shells) with 
equal probability. This simple model is favored for its 
ability to produce a relatively high radiation efficiency 
and a narrow peak energy distribution within a same 
burst. The new ingredient we consider here includes the 
ri-component for both the fast and the slow shells. 

Comparing the prompt emission energy and the after- 
glow kinetic energy derived from multi-wavelength data 

^ In this work, we only discuss the purely hydrodynamic fireball 
model. If the outflow is magnetohydrodynamic, the n-p decoupling 
process is different (Vlahakis et al. 2003). 
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fits (Panaitescu & Kumar 2002), it is found that a signifi- 
cant fraction of the initial kinetic energy is converted into 
internal energy and radiated as gamma-rays (Panaitescu 
& Kumar 2002; Lloyd-Ronning & Zhang 2004). Within 
the internal shock model, this requires that the velocity 
difference between the shells is significant, i.e. ryf 3> ?7s, 
and that the two masses are comparable, i.e. mt « rrig 
(Piran 1999). Generally, rjs is in the order of tens, and 
r]{ is in the order of hundreds. Thus, for slow shells 
the n, p components likely coast with a same LF, i.e. 
Ps,!! — Ps.p — Vsj while for fast shells the n, p compo- 
nents may have different LFs, which are denoted as Ft^n, 
Pf^p, respectively. 

When an inner faster shell catches up with an outer 
slower shell, the ion component merges into a sin- 
gle ion shell. For simplicity this process is approxi- 
mated as a relativistic inelastic collision. Energy and 
momentum conservations result in Fm — [(wf^pPf^p -|- 

Ws,pPs,p)/(TOf,p/Pf,p -I- ms.p/Ps.p)]^/^ (Paczynski & Xu 
1994; Piran 1999), where P,„ is the LF of the merged ion 
shell. 

As shown in Fan & Wei (2004), the merged ion shells 
are further decelerated by the decay products of the slow 
neutron shells. Given a small LF of the slow neutron 
shells, the mean neutron decay radius is smaller, i.e. ~ 
900cPs,„ ~ 8 X lO^^'Fi.s cm. The decay products would 
be collected by the merged ion shells at the radii in the 
range oi 10^^ — lO^^cm. Shocks are formed, and the 
ion shells further merge with the decay products of the 
slow neutron shell. The resulted typical LF of the ion 
shells would be 



^ V ('"f:P + "^s,p)/Pm + ms,n/Ps,n ' 

which will be regarded as the initial LF of the ion shell 
(hereafter called the I-ejecta) whose dynamical evolution 
we will discuss below. The thermal LF of the protons in 
the shocked region reads 

V (™f,P + '^s,p)rm + TOs.nPs.n 

7th,i ~ : : 

mt,p + ms,p + ms,n 

\f (TOf,p -I- ms,p)/Pm + ms,n/Ps,n - 1 , (2) 

and the electron synchrotron emission in the region gives 
rise to a bright UV flash (Fan & Wei 2004). 

Notice that the presence of the neutron component 
does not help to increase the 7-ray emission efficiency. 
Given a same amount of baryon loading and a same to- 
tal energy, the LF of the fast ion shells is larger than that 
in the neutron-free model, while the LF of the slow ion 
shells (~ tens) remains the same. Although the collision 
of the fast-slow ion shells is more efficient than in the 
neutron free model, a significant part of the total energy 
is carried by the fast neutrons shell (i.e. the N-ejecta) 
most of which is not translated into the prompt 7— ray 
emission. As a result, the GRB efficiency in the neutron- 
fed model may be even lower than that in the neutron free 
model. Rossi et al. (2004) gave another argument on the 
low radiation efficiency of neutron-rich internal shocks. 
It is worth noticing that the internal shock efficiency is 
also lowered in a magnetized flow, even magnetic dissi- 
pation plays an essential role (Fan et al. 2004c). 



2.3. The trail of the N-ejecta 

While the fast ion shells are decelerated by the slower 
ions (including both the slow ion shells and the decay 
products of the slow neutrons), the fast neutron shells 
(which will be denoted as the "N-ejecta") are not. As 
long as they move fast enough, say, Fn = Pf,n > Pi 
(Pi is the LF of the I-ejecta, considering the dynamical 
evolution), the N-cjecta would freely penetrate through 
all the ion shells in front of them and would separate from 
the I-ejecta more and more. As the result of /3-decay, the 
mass of the fast neutrons gradually decreases as 

Af„(i?) = M^eM-R/Rp), (3) 

where M° is the original mass of the total fast neutrons, 
Rp = PnCT is the characteristic neutron decay radius, 
and T ~ 900 s is the rest-frame neutron decay time 
scale. If neutrons decay in the circumburst medium, the 
/3— decay products p and e~ would immediately share 
their momenta with the ambient particles through two- 
stream instability (B03). Here we take the number den- 
sity of the external medium (in unit of cm~^) as 

r const., (ISM), 
" ^ \ 3.0 X lO^^A.R-^, (wind), 

where "ISM" and "wind" represent the uniform inter- 
stellar medium and stellar wind, respectively, A* = 
(M/10-^Mq yr-i)(z;w/103km s-i)-\ M is the mass 
loss rate of the progenitor, and Ww is the wind veloc- 
ity (Dai & Lu 1998; Meszaros et al. 1998; Chevaher & 
Li 2000). 

B03 discussed a neutron decay "trail" , which is a mix- 
ture of the ambient particles with the decay products left 
behind the neutron front. He then discussed the interac- 
tion of the follow-up I-ejecta with this trail and suggested 
that the trail would modify the dynamical evolution of 
the I-ejecta and hence the afterglow lightcurves. Such 
a clean picture is strictly valid for very short bursts, 
for which the N-I ejecta separation radius Rs fea. j34|'l 
is shorter than R/j. For most of the long GRBs for which 
the majority of the afterglow data are available, not all 
the decayed mass |dMn| is deposited onto the circum- 
burst medium. Part of it is deposited onto the I-ejecta 
itself. The neutron decay products therefore influence 
the dynamical evolutions of both the medium and the 
ejecta. In this paper, we will define the "neutron trail" 
only as those decay products that are deposited onto the 
circumburst medium. Those deposited onto the I-ejecta 
itself will be treated separately. 

2.4. The dynamical evolution of the system 

With the presence of neutrons, the dynamics of the 
whole system is far more complex than the neutron-free 
case. The physical pictures for both the ISM and the 
wind environment are considerably different. Below is a 
qualitative discussion, and more detailed formulation is 
presented in f|3| 

For the wind case, the medium is very dense, so the 
decay products are immediately decelerated significantly, 
and the resulting bulk LF of the trail is in the order of 
tens. The trail is picked up by the I-ejecta quickly. A PS 
propagates into the trail and a RS propagates into the 
I-ejecta. Since the LF of the FS F^g is smaller than Pn, 
the FS never propagates into the wind medium directly. 
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although at later times the neutron decay rate becomes 
so low that the neutron trail essentially does not alter 
the wind medium. Besides the emission from both the 
FS and the RS, for x '^-^ 1 the unshocked neutron trail 
may also give non-negligible emission contribution in the 
optical band. The thermal LF in the trail is several for 
typical parameters. Given a reasonable estimate of the 
local magnetic fields, the typical synchrotron radiation 
frequency is around the optical band, so that it could 
give an interesting contribution to the optical afterglow 
emission. This is discussed in H3.1l in detail. 

For the ISM case and for i? < a few Rp, unless n is 
very large, the inertia of the swept medium is too small 
to decelerate the decayed fast neutrons significantly. If 
the I-ejecta moves slower than the N-ejecta (which is 
our nominal case), the neutron trail most likely moves 
faster than the I-ejecta, i.e. 7 > Fj, where 7 is the LF 
of the trail. A gap forms between the trail and the I- 
ejecta. Since the neutron decay rate drops with radius, 
the trail LF also decreases with radius. This leads to 
pile-up of the neutron trail materials with higher LFs. A 
self-consistent treatment of the dynamical evolution of 
the trail is complicated. For the purpose of this paper, 
it is adequate to treat the fast trail as a global T-ejecta, 
which in many aspects are analogous to the I-ejecta. In 
front of the T-ejecta, there are still some newly decayed 
products (we still call these the neutron trail) behind 
the N-ejecta front. Even further ahead the N-ejecta is 
the unperturbed ISM. The fast moving T-ejecta shocks 
into the trail and further into the ISM and gets decel- 
erated. The Lorcntz factor of this forward shock F^g is 
sometimes larger than the Lorentz factor of the N-ejecta 
Fn, so that the shock front is directly in the ISM. In 
any case, the very early afterglow emission is from this 
"T-ejecta" forward shock. The I-ejecta, which lags be- 
hind without significant interaction with the ISM, is not 
noticeable in the beginning. Later, the I-ejecta catches 
up with the decelerated T-ejecta, giving rise to a pair 
of refreshed shocks that power strong IR/Optical emis- 
sion. An energy-injection bump is expected, which has 
direct observational consequence. The whole process is 
discussed in i j3.2l in detail. 

3. EARLY OPTICAL AFTERGLOW LIGHTCURVES: LONG 

GRBS 

Below we calculate the R-band early afterglow 
lightcurves, focusing on the novel properties introduced 
by the neutron component. In this section we deal with 
the traditional long GRBs, those with T90 > 2s and be- 
lieved to be produced during core collapses of massive 
stars. Short bursts are discussed in 21 Two widely 
discussed circuniburst medium types are ISM and wind, 
and we discuss the lightcurves for both cases in turn. 
It is widely believed that synchrotron radiation of the 
shocked electrons powers the observed afterglow emis- 
sion. The detailed formulae for synchrotron radiation 
are summarized in Appendix ^ In the text we mainly 
focus on the dynamical evolution of the entire system. 

One novel feature of the neutron-fed early afterglow is 
the emission from the neutron decay trail, especially in 
the wind case. The detailed particle acceleration mech- 
anism as well as the radiation mechanism in the trail 
is poorly known. In this paper, we adopt two extreme 
models to treat the trail emission (see Appendix for 



details). The first treatment is similar to the shock 
case, i.e., assuming a power-law distribution of the elec- 
trons. The second treatment is to simply assume a mono- 
energetic distribution of the electron energy. We believe 
that a more realistic situation lies in between these two 
cases. In both cases we introduce the equipartition pa- 
rameters of electrons and magnetic fields. 

Similar to Fan & Wei (2004), we take Etot = 2.0 x 
2^Q53gj.gg (including the total energy of the fast/slow 
neutrons and protons at the end of the prompt 7— ray 
emission phase), Fn = 300, Fg^n = 30, Fm = 200, 
'7if,p = ™s,p- In the ISM and the wind cases, the num- 
ber density is taken as n = 1 and n = lO^^i?"^ (i.e. 
A* = 1/3), respectively. The width of the N-ejecta 
in the observer frame is taken as A = lO^^cm, where 
A — cTgo/{l + z) is the width of the ion-ejecta in the 
rest frame of the central engine and the observer (Cor- 
respondingly, T90 = 33(1 -I- 2)s). In order to find out 
how the neutron emission signature depends on the neu- 
tron fraction, in each model we calculate four lightcurves 
that correspond to the neutron-to-proton number ratio, 
X = 0.0, 0.1, 0.5, 1.0, respectively. 

Besides the neutron signature, another important in- 
gredient in the early afterglow phase is the RS emission. 
The early dynamical evolution of a neutron- free shell and 
its interaction with the circumburst medium has been in- 
vestigated in great detail both analytically (e.g.. Sari & 
Piran 1995; Kobayashi 2000; Chevalier & Li 2000; Wu et 
al. 2003; Z. Li et al. 2003; Kobayashi & Zhang 2003b; 
Fan et al. 2004a; Zhang & Kobayashi 2004) and numer- 
ically (Kobayashi & Sari 2000; Fan et al. 2004b; Nakar 
& Piran 2004; Yan & Wei 2005; Zou et al. 2005). The 
early afterglows of several GRBs (including GRB 990123, 
GRB 021211 and possibly GRB 021004) have been mod- 
eled within the RS emission model (e.g.. Sari & Piran 
1999; Meszaros & Rees 1999; Wang et al. 2000; Soder- 
berg & Ramirez-Ruiz 2002; Fan et al. 2002; Kabayashi 
& Zhang 2003a; Wei 2003; Zhang et al. 2003; Kumar & 
Panaitescu 2003; Panaitescu & Kumar 2004; McMahon 
et al. 2004). With the presence of neutrons, the dynam- 
ical evolution of the ejecta and the RS emission becomes 
more complicated. Below we will formulate the entire 
process in detail. The results are reduced to the conven- 
tional results when x = 0, but become more complicated 
as X gets larger. 

In this work we assume that the medium contains elec- 
trons and protons only. In reality, at radii smaller than 
lO^^cm, the medium may be enriched by pairs cre- 
ated by the interaction of the prompt 7— rays with the 
back-scattered 7-rays by the medium (e.g. Beloborodov 
2002 and references therein). The influence of pairs on 
the early afterglow lightcurves is neglected in this work 
and will be considered elsewhere. 

3.1. Wind case 

We will first consider the wind model in which no gap 
is formed between the I-ejecta and the trail so that the 
dynamics is simpler (see ^3. 21 for a detailed treatment for 
the gap case). The ejecta-trail interaction can be divided 
into two phases, (i) Before the RS crosses the ejecta, 
R < Rx (where Rx is the radius at which RS crosses 
the ejecta), there exist two shocks, i.e. a FS expanding 
into the trail, and a RS penetrating into the I-ejecta. (ii) 
After the RS crosses the ejecta, R > R^, only the FS 
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Fig. 1. — A schematic diagram of a neutron- fed long GRB in the 
wind interaction case. 



exists. 

3.1.1. The dynamical evolution for R < Ry 

As shown in Figure^ there are five regions: (1) the rest 
wind medium; (2) the unshocked neutron trail, moving 
with a LF 7; (3) the shocked neutron trail, moving with 
a LF F; (4) the shocked I-ejecta, moving with the LF 
F; and (5) the unshocked Fejecta, moving with a LF Fj. 
Notice we have generally defined the last LF as Fj rather 
than Fj since the decay products directly deposited onto 
the I-ejecta would change Fj from its initial value. 

We first treat the dynamics of the neutron trail. Fhe 
main purpose is to calculate the velocity of the mixture of 
the ambient particles with decay products (i.e. the trail) 
P{R) = v{R)/c. Following B03, we take dm, \dMn\ as the 
mass of the ambient medium overtaken by the N-ejecta 
and the mass of the decayed fast neutrons, respectively. 
For the long bursts we are discussing, only a fraction of 
\dMn\ is deposited onto the medium (the rest is deposited 
onto the Fejecta). A shock propagates into the trail, 
and in the problem what is relevant is the fraction of the 
decay products that are deposited onto the unshocked 
trail (region 2 in Figure QJ. We denote this fraction as 



/ = mm{l, — - — }, 



(5) 



where A is the width of the N-ejecta, which is the same as 
that of the Fejecta, and An-f is the separation between 
the front of the N-ejecta and the front of the FS that 
propagates into the neutron decay trail. This is valid for 
the wind case, in which An-f > is always satisfied. 
For the ISM case f tl3.2|l . under certain conditions the 
FS front leads the N-ejecta front, and we take An-f < 
following the same definition here. If this happens, 
nothing is deposited onto the unshocked medium, and in 
our treatment we will define f — whenever An-f < 
is satisfied. Notice that by defining /, we have assumed 
that neutrons in the N-ejecta decay uniformly within the 
length A, which is justified by the fact A R. The same 
applies to the other two fraction parameters (g, ea.|16j 
and h, eq.021 defined below). 



The width An-f is determined by^ 

dAN_F = (/?r„ -/3rps)rf^, (6) 
where Fps is the LF of the FS and hereafter we denote 
/3rA a^s the dimensionless velocity corresponding to the 
LF Fa. 

Considering an inelastic collision between the medium 
(with mass dm) and the decay products (with mass 
f\dMn\), energy and momentum conservations give 

/FnMMnI + dm = 7(7th + l)(dm + f\dAQ), 

/F„/3rJrfMn|=7/37(7th + l)(fim + /|(iM„|) , (7) 

where 7th is the thermal LF of the mixture (exclude the 
rest mass). With equations iQ, we have (see also BOS) 



/3^(i?) 



1 + (FnC)- 



7(i?) 



(C2 + 2F„C+ 1)1/2' 



where we have defined a parameter 
.\dAUR)\ 



dm{R) 



Mn dm 
■' RpdR' 



fMl , R 



(-f )i-^exp( 



R_ 

Rji ' 



(8) 



(9) 



which denotes the deposited neutron mass on unit 
medium mass (notice that it is slightly different from 
the same parameter defined in BOB with the correction 
introduced by the / factor). Here is the mass of the 
medium contained within R < R^, and k is a parame- 
ter denoting the type of the medium, i.e. k = I for the 
wind case and fc = 3 for the ISM case. For the wind case 
discussed here, we have 



m)] 



■exp( 



R_, 

Rp ' 



(10) 



Equations lO also yields the initial thermal LF of the 
neutron trail 

7th = (C' + 2F„C + 1)'/V(1 + 0-1- (11) 
For Fj^i <C C Tn (which is valid for the wind case), 
the internal energy of the trail far exceeds its rest mass 
energy, i.e., 7th ^ 1, so that significant radiation is ex- 
pected. When C, > Fn (e.g. the ISM case), one has 
7th <C 1 so that the trail emission is unimportant. 

The total proton number density (in its proper frame) 
in the unshocked trail, including the neutron decay prod- 
ucts, can be expressed as (B03) 

ntr = n(l-f C)(C' + 2F„C+1)'/'. (12) 
We define Ap-R as the width of the shocked regions, 
i.e. the distance between the FS and the RS. One then 
has 

dAp^n = iPr^s - 0r^s)dR, (13) 

where /Spps /^rj^g are the dimensionless velocities of 
the FS and the RS, respectively, which read (e.g.. Sari & 
Piran 1995; Fan et al. 2004b) 

r/3r(4723 + 3)-7/37 



Trs ' 



r(4723 + 3) - 7 ' 
r/3r(4745 + 3) - Fi/3r, 
F(4745 + 3) - Fi 



(14) 



"* Hereafter we have adopted dR ~ cdi, where i is the time in 
the rest frame of the emission source, so that the time interval is 
replaced by dR. 
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where 723 is the LF of the region 2 relative to the region 

3, and 745 is the LF of the region 5 relative to the region 

4. They are calculated as 

723 = 7r(l - Pjl3r), 745 = rir(l - Pr^Pr)- (15) 

We define another fractional parameter g as the frac- 
tion of the decay products that are deposited onto the 
shocked region (both regions 3 and 4 in Figure^. We 
then have 
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Af-r/A, for An-r < A 
1 - /, for An-r > A 



(16) 



where An-r is the separation between the front of the 
N-ejecta and the RS front, which is described by 



dAN-R = (/3r„ -/3rRs)rfi?- 



(17) 



When An-r < A the RS has not passed the end of the 
N-ejecta, while when An-r > A, the RS passes through 
the end of the N-ejecta, reaching the region where no 
neutron decay product is being deposited. 

For An-r < A, some decay products are being de- 
posited onto the I-ejecta directly (region 5 in Fig^, 
which will alter the dynamics of the I-ejecta. We can 
treat the process in the way similar to the treatment 
of the neutron trail onto the circumburst medium. The 
difference here is that the I-ejecta itself is highly rela- 
tivistic. The energy and momentum conservations are 
more analogous to those to describe internal shock colli- 
sions. Unlike the circumburst medium case in which the 
trail LF only depends on R, one needs two parameters 
to describe the I-ejecta dynamics. Besides the radius R, 
one needs to specify the location of the I-ejecta element 
in the I-ejecta proper. We parameterize this by the dis- 
tance from the initial I-ejecta front, which is denoted as 
Ai so that < Ai < A is satisfied. After some algebra, 
the LF of the I-ejecta layer with the coordinate Ai at a 
particular R could be written as 



ri(Ai,i?) 



Cr„ + (i + 7tVi)r? 



^C7r„ + (i + 7li)/rr 

and the corresponding thermal LF reads 



(18) 



7th,i(Ai,i?) 



'(Cr„ + r?)[c/rn + (i + 7tVi)/r?] 



where 



C = C(Ai,i?) 



(1 + C)^ 



1 — exp 



' Rp 



(19) 



(20) 



and is the initial mass of the I-ejecta. The parameter 
C, is analogous to the C parameter (eq.[3) for the trail 
treatment, which denotes the average decayed neutron 
mass on unit I-ejecta mass. The characteristic radius 
Ri depends on whether or not the I-ejecta layer (with 
coordinate Ai) overlaps the N-ejecta as the reverse shock 
sweeps the layer (so that the layer no longer belongs to 
the region 5). If the layer still overlaps with the N-ejecta 
as it is swept by the RS, through out the period when 
the layer is in the region 5 there has been always neutron 
decay products deposited onto the layer. In such a case, 
one simply has i?i = R. For those layers that the end 



of the N-ejecta has passed over earlier so that there is 
no neutron decay deposition at the RS sweeping time, 
one has to record the decayed neutron amount at the 
final moment with N-ejecta overlap. In this case, one 

has i?i « 2(r?r„)2(A - Ai)/(r2 - vf). 

The RS crossing radius i?x can be calculated according 

to 

dAi-R = A, (21) 

/o 

where Aj-r is the distance between the I-ejecta front and 
the RS assuming no shock compression. It is determined 

by 

dAi-R = (/3r: -/3rKs)rf^- (22) 

With the above preparation, one can delineate the dy- 
namical evolution of the system for R < i?x using the 
following equations. After some algebra, the energy con- 
servation of the system yields 

(M + U/c')dr = (7 - r723)(l + 7th)rfmt,. + (Fi 
-r745)(l + 0(1 + 7th,i)dmi + 5(r„ - r)|dA/„|, (23) 

in deriving which, the small thermal energy generated 
by the deposition of the neutron decay products into the 
shocked region has been neglected. Here U is the thermal 
energy generated in the two shock fronts, which can be 
calculated by 

dU/c^ - (1 - ei)[723(l + 7th) - l]dmtr 

+ {1- e2)(l + C)[745(l + 7th,i) - l]rfmi, (24) 

where ei and £2 are the radiation efficiency (see the Ap- 
pendix ^ for definition). The differential increase of the 
trail mass, the I-ejecta mass (just initially) and the total 
mass in the shock region can be expressed as 



dmtr = [1 + C{R)]ATmR'^dR, 
dmi = Mi°(iAi-R/A, 
dM = dmtr + {1 + C)dmi + g\dMn\ 



(25) 
(26) 
(27) 



Combining equations (|13I27|) with the well known 
radius-time relation 



dR = ^FpsC dt/{l - Ptps)^ 



(28) 



one can calculate the dynamical evolution for R < Rx 
(see our numerical example in H3.1.4|l . Hereafter t — 
iobs/(l + z) denotes the observer's time corrected for the 
cosmological time dilation. 

3.1.2. The dynamical evolution for R > i?x 

After the RS crosses the I-ejecta, only the FS exists. 
Equations (|23l 1241 127|l can be simplified as 

{M+U/c')dT = (7-r723)(l+7th)dmt,+<?(r„-r)d|M„|, 

(29) 

dU/c^ - (1 - ei)[723(l + 7th) - l]rfTOtr, (30) 
dM = dnitr + g\dMr,\. (31) 



7 




Fig. 2. — The early optical afterglow lightcurves of a neutron- 
fed long GRB in the wind interaction case, (a) The dynamical 
evolution of the LF of the shocked region as a function of time, (b) 
R-band lightcurves, with the IC cooling effect due to the prompt 
7— rays interacting with the shocked regions being ignored. Thick 
lines include contributions from all emission components, includ- 
ing the FS, RS and the neutron decay trail. Thin lines are for 
trail emission only (for a power-law energy distribution of the 
electrons, see Appendix |^ for details). The dotted, dash-dotted, 
dashed and solid lines represent x = 0.0, 0.1, 0.5, 1.0 respec- 
tively. Following input parameters are adopted in the calculations: 
Etot = 2.0 X 10"^ ergs, A = lO^^cm, z = l [i.e. Dl = 2.2 X 10^* cm, 
which corresponds to the standard (Qrn,^A) = (0.3,0.7) ACDM 
cosmological model], Fn = 300, Fm = 200, Fs.n = 30, and 



lO^Sc 



(i.e. j4t = 1/3), respectively. The parame- 



ters Ce = 0.1, EB = 0.01 and p = 2.3 are adopted for the FS and 
RS shocks as well as the trail, (c) R-band lightcurves, but the 
IC cooling effect due to the prompt 7— rays overlapping with the 
shocked region and the trail has been taken into account. The av- 
eraged 7— ray luminosity is taken as = 1051 ergs s 1 . Other 
parameters and line styles are the same as those in (b). 



3.1.3. The impact of the initial prompt j—rays on 
electron cooling 

As pointed out by Beloborodov (2005), if the RS emis- 
sion overlaps with the initial prompt 7— ray emission, 
the cooling of the shocked electrons is likely dominated 
by the inverse Compton (IC) scattering off the initial 
prompt 7— rays. This is a very common case for the wind- 
interaction case, as has been studied in detail in Fan et 
al. (2005a), who discuss the interesting sub-GeV photon 
emission and high energy neutrino emission due to the 
overlapping effect. In the ISM case, usually the prompt 
7— ray emission has crossed the RS region before the RS 
emission becomes important, so that the above overlap- 
ping effect is not as important as in the wind model^. 
We have performed calculations of the overlapping effect 
in the ISM model as well, and found that the resulted 
early optical afterglow lightcurves are nearly unchanged 
for typical parameters. So we do not present the result in 
FigurelHl Nonetheless, we would like to point out that for 

5 We can understand this as follows. For illustration, we take 
the neutron-free model (x = 0). The ejecta are decelerated signifi- 



1/3 



-1/3 



r'^/ t , at which 



cantly at the radius -Rdcc ^ lO^^^cm 
the rear of the prompt 7— ray emission leads the rear of the ejecta 
at an interval ~ i?dec/2r^ ~ 1.2 X lO^^cmEV^g ^''^F;;*2'' 3, 
which is usually larger than (or at least comparable to) the width 
of the fireball A ~ lO^^^cm. As a result, for the typical parameters 
taken in this work, the IC cooling due to the overlapping in the 
ISM case is not as important as in the wind model. 



some GRBs born in the ISM environment whose duration 
is long enough and/or whose bulk LF is large enough, the 
overlapping effect should be also very important. This 
has been the case of GRB 990123 (Beloborodov 2005), 
which has a long (Tgo ~ 120s) duration and a large LF 
(- 1000). 

Similar to Fan et al. (2005a), we assume that the in- 
ternal shock efficiency is r ~ 0.2 and define the IC cool- 
ing parameter Ys = U^/Ub- Here U-y ~ L-^/AnR'^r'^c is 
the MeV photon energy density in the rest frame of the 



shocked region, and Ub — B 



{TyTfeBTiAmpC^ is 



the co-moving magnetic energy density in the same re- 
gion, where ~ rci<^tot/(l — f)^ is the luminosity of 
the 7— ray burst, and 714 « Etot / A-k R^T^^ m-pC^ IS. is the 
co-moving density of the un-shocked ejecta. After some 
simple algebra, we derive Ys ~ ^/[{l — r)eB]. Notice that 
this is only valid when the Klein-Nishina correction is 
unimportant, i.e. a; <C 1, where x = e^jc/^n^cC^ 7c/r, 
e-y ~ TOqC^ is the energy of the prompt 7— rays, jc is the 
random LF of the emitting electrons. More generally, 
we should make the Klein-Nishina cross section correc- 
tion, i.e. aic = A{x)aT, where A(x) = |[^{^^^ - 

"^^^^ ^^'^ asymp- 
2x + 2^ for X < 1, and 



ln(l + 2x)} + ^ln(l + 2x) 
totic limits A{x) w 1 



Aix) 



-^{\n2x i) for a; > 1 (e.g. Rybicki & Light- 



man 1979). Notice that such a correction is necessary 
since the target photons, i.e. the initial prompt 7— rays 
are usually energetic enough in the rest frame of the elec- 
trons. As a result, for the electrons accelerated by the 
reverse shock, for x ^ 1, the Compton parameter reads 
(see Footnote 1 of Fan et al. [2005a] for derivation) 



A{x)r 
x{l - r)eB 



(32) 



Usually, the magnetic field densities in the FS and the 
RS are nearly the same (e.g. Sari & Piran 1995), so Eq. 
H32I) also applies to the FS shock region. The overlapping 
of the initial prompt 7— rays with the FS emission lasts 
until / (1 - l3Ys)dR = A. 

Similarly, for the electrons accelerated in the neutron 
trail, we can also introduced the IC cooling parameter 
Yt ~ A{xi)U^^t/ {xiUb'^^) for xi ^ 1, where U^^t ~ 
L^l A'kR^^'^c, is calculated by Eq. (jBip . and xi = 
e^-fchm^c^ ~ 70/7. 

When the prompt 7— rays overlap the shocked regions 
and the trail, to calculate the coolin g fre quency, the Y 
parameter in Eqs. ljA3|l . (jB3|) and l|B4p should be re- 
placed by Ys (for electrons accelerated from the shock) 
or It (for electrons accelerated in the neutron trail). 

3.1.4. Model lightcurves 

In Figure |21 we present the dynamical evolution of 
the system (Figure El^a)) and the early optical after- 
glow lightcurves in Figs|2Ib) and (c). In Figl^^b) the 
IC cooling effect due to the initial prompt 7— ray over- 
lapping with the shocked region and the trail has been 
ignored. In FigI2tc), this IC cooling contribution is ex- 
plicitly taken into account. The neutron-to-proton num- 
ber density ratio x is adopted as x = 0.0, 0.1, 0.5 and 1, 
respectively. 

As shown in Figl^^a), for i? < i?x , T is nearly a con- 
stant. This is consistent with the familiar result in the 
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neutron- free wind model (e.g., Chevalier & Li 2000; Wu 
et al. 2003; Kobayashi & Zhang 2003b; Fan et al. 2004b; 
Yan & Wei 2005; Zou et al. 2005). The main reason 
is that the density of the N-ejecta also drops as cx R~^, 
and that the fraction of the decayed neutrons essentially 
remains the same during the shock crossing process. The 
latter effect is manifested by the fact that Rx is compara- 
ble to Ri3 for the current nominal parameters, so that the 
neutron decay rate [oc exp(— remains essentially 
the same for R < i?x • 

The dynamical evolution of the shocks presented here 
[e.g. Fig|2Ia)] look different from the one presented 
in Fig.l of B03, which shows that the dynamics of a 
neutron-fed fireball is very different from the neutron- 
free one. This is because B03 presented the dynamical 
evolution of the trail in the radiative regime, i.e., the 
total thermal energy of the trail is radiated promptly be- 
fore the trail is picked up by the I-ejecta. Here we assume 
an electron equipartition parameter Cc = 0.1 in the neu- 
tron trail^. For such a small Eo, even in the fast cooling 
regime, only 0.1 of the total thermal energy is radiated, 
and the system could be still approximated as an adia- 
batic one. Our results indicate that as long as the trail 
is in the quasi-adiabatic regime, the bulk of the energy 
is still dissipated in shocks, and the existence of the neu- 
trons does not influence the dynamical evolution of the 
system significantly (see also BOS for similar discussions) . 

In any case, we include the emission contribution from 
the neutron trail by assuming a power-law distribution 
of the electrons, and find that it gives an interesting sig- 
nature when X is close to unity. In Fig[21[b) we present 
the lightcurves for different x values, with the IC cool- 
ing by the initial prompt 7— rays ignored. The thick 
lines are the synthesized lightcurves including the con- 
tributions from the FS, RS and the neutron trail. In 
order to characterize the function of the trail, we plot 
the trail emission separately as thin lines. We can see 
that the early afterglow for x — ^ shows a bright plateau 
lasting for ~ 100s, and the main contribution is from 
the trail. For x ^ 0.5, an early afterglow plateau with 
a shorter duration is still evident, but it is mainly due 
to the contribution of the RS emission^. Comparing the 
lightcurves with small x's (e.g. 0.0 and 0.1) to those with 
large x's (e.g. 0.5 and 1.0), one can see that the early 
afterglow intensity is much stronger for high-x case, al- 
though the dynamical evolution of the I-ejecta is rather 
similar [FigI21^a)]. The reason is that the trail deposits 
an electron number density about 10 times of that in the 
wind medium, so that the total number of the emitting 
electrons in the shocked region is greatly increased. 

In Fig[2Ic) we present the lightcurves by taking into 

® As noted in BOS, the energy dissipation mechanism in the neu- 
tron decay trail could be fundamentally different from that in the 
coUisionless shock. Nonetheless, for known mechanisms of particle 
acceleration in a proton-electron plasma, the fraction of thermal 
energy given to electrons is usually significantly smaller than that 
given to protons. We therefore suspect that Ec ~ 0.1 may be rea- 
sonable. 

''' The lightcurve temporal index before shock crossing is ~ 
rather than 1/2 discussed in other papers (e.g. Chevalier & Li 
2000; Wu et al. 2003; Kobayashi & Zhang 2003b; Fan et al. 2004b; 
Yan & Wei 2005; Zou et al. 2005). The main reason is that for the 
nominal parameters in this paper, the optical band is above both 
i^m and Uc for the RS emission, Fi,^ oc t~(P~2)/2^ which is very 
flat. 



account the IC cooling effect due to the initial prompt 
7— rays overlapping with the shocked regions and the 
trail. Line styles are the same as in Fig. ^h). It is 
apparent that the very early R-band emission, including 
the FS emission, the RS emission and the trail emission, 
has been suppressed significantly. This means that the 
IC cooling effect is very essential in shaping the early 
afterglow lightcurves in the wind case. After the rear 
of initial prompt 7— ray front has crossed the FS front, 
the FS emission becomes rather similar to that of the 
FigEJb). 

In the wind case, the difference between the neutron- 
rich lightcurve and a neutron-free lightcurve is only 
marginal (both Fig[2Ib) and Fig|2Ic)). One potential 
tool to diagnose the existence of the neutron component 
is to search for the trail emission component. This is in 
principle not straightforward since there are many un- 
certainties in categorizing the trail emission (Appendix 
p. In Fi gure 13 (with the IC cooling effect ignored), 
we calculate the expected trail emission for two extreme 
models for electron energy distribution, i.e. a power law 
model and a mono-energetic model. A more realistic 
model should lie in between these two models. Figure 
O indicates a fortunate fact that the differences between 
the two extreme models are far from large, although the 
mono-energetic model predicts a stronger R-band trail 
emission. This gives us confidence that the crude treat- 
ment of the trail emission in Fig|21 gives a first-order pre- 
sentation of the reality. Another remark is that differ- 
ent electron energy distributions result in different trail 
emission spectrum. In particular, a mono-energetic or 
thermal-like distribution gives rise to a much lower emis- 
sivity at high frequencies. Even in the optical band, 
multi-color observations at early times can be used to 
diagnose the existence of the neutron trail. For example, 
if the trail electron energy distribution is thermal-like or 
in a form largely deviated from the simple power-law, 
multi-color observations around 100s with the UVOT on 
board Swift (there are 6 colors in the 170-650 nm band 
for UVOT) can reveal interesting clues of the existence 
of neutrons in the fireball. 

3.2. ISM case 

We now turn to discuss the ISM-interaction for long 
GRBs. The main characteristic of the ISM case is that 
the trail typically moves faster than the I-ejecta, so that 
a gap forms between the I-ejecta and the T-ejecta (the 
trail ejecta). The latter interacts with the ISM and dom- 
inates the early afterglow. Later the I-ejecta catches up 
and produces an energy injection bump signature. Below 
we will calculate the typical lightcurves with the nomi- 
nal parameters adopted in this paper. Two stages will 
be considered separately, i.e. (1) when the T-ejecta is 
separated from the I-ejecta and dominates the afterglow 
emission, and (2) when the I-ejecta collide onto the decel- 
erated T-ejecta. Except otherwise stated, the notations 
adopted in this section are the same in H3.ll when the 
wind case is discussed. 

3.2.1. The dynamical evolution of T-ejecta before 
I-ejecta collision 

The C{R) parameter (eq.PI) in the ISM case reads 

[C(i?)pM^Zi^(^)-2exp(-A), (33) 
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Fig. 3. — The trail emission lightcurve in the wind interac- 
tion case, with the IC cooUng effect due to the initial prompt 
7— rays overlapping with the trail ignored. The dotted, dashed 
and solid lines represent x = O-l; 0-5, 1-0 respectively. The thick 
lines (marked by P) represent the emission by assuming a power 
law distribution of electrons, and thin lines (marked by M) repre- 
sent the emission by assuming a mono-energetic distribution. The 
initial input parameters are the same as those presented in Figure 
I^Jb) caption. 
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For X = li the nominal parameters are Fj — 100, 
Fn = 300, -f°^j = 0.433 and CiR) = 2.85 x 
W^fiR/Rpy^expi-R/R/^). The N-ejecta leads the F 
ejecta, and completely separates from the Fejecta at 



Rs 



A 



- Pro (F2 



(34) 



As long as f exp{-R/ Ri3){R/ R^)-^ > 2.34 x lO^^, one 
has 7(i?) > Fj (eq.|S]). The Fejecta lags behind the trail, 
and a gap forms between the Fejecta and the trail. This 
happens for x > several x 0.1, when the decay products 
deposited onto the ISM have a large enough momentum 
to drag the trail fast enough to lead the Fejecta. With 
equation H34I) , we can estimate how many decay products 
have been deposited into the Fejecta directly 

:M„^ / (1 



R ^ ^ i? ^ dR 

Rs Rp Rf3 



= [1 



exp(— a) 1 



(35) 



where a = Rs/Rp- For a > 2, we have M^^d > 0.5Af°, 
i.e., most of decay products have been deposited into the 
Fejecta. 

For the trail, the layers from the behind moves faster 
than the leading layers, since [C(^)]^^'^ decreases with 
R. This leads to pile-up of the trail materials. The faster 
trail materials would shock into the slower trail materials 
in front. In the following treatment, we approximately 
divide the trail materials into two parts. The fast trail 
(including forward shocked trail) is denoted as the T- 
ejecta, while the upstream unshocked part is still called 
the trail. The location of the shock is roughly defined 
by requiring the bulk FF of the T-ejecta (which is cal- 
culated through energy and momentum conservation) is 
larger than 27, so that the relative FF between the two 
parts exceeds 1.25. Since the trail is hot with a relativis- 
tic temperature, the upstream local sound speed c/\/3 



Fig. 4. — The schematic diagram of a neutron-fed long GRB 
in the ISM-interaction case, (a) When the T-ejecta (region 3) is 
separated from the I-ejecta and interacts with the trail or ISM. (b) 
When the l-ejecta collides into the T-ejecta producing refreshed 
shocks. 



corresponds to a local sound FF ~ •\/3/2 ^ 1.22. So 
our shock forming condition is self-consistent. The shock 
front moves faster than the T-ejecta, and Ptfs 
ceed /3r„- If this persists long enough, the N-ejecta lags 
behind the FS front. The decay products deposit into 
the shocked material directly rather than mixing with 
the ISM. The shock propagates into the ISM directly. In 
such a case, the dynamics is simplified. 

Before collision, there are six regions in the whole sys- 
tem (see FiglUa) for illustration): (1) the rest ISM; (2) 
the unshocked trail moving with 7; (3) the T-ejecta (in- 
cluding the shocked trail between the FS and the RS) 
moving with F; (4) the gap between the T-ejecta rear 
and the l-ejecta front, in which the decayed protons move 
with F„; (5) the overlapped region of the l-ejecta and 
the N-ejecta; and (6) free moving l-ejecta. 

The width of the T-ejecta (i.e., region 3), Af-r is still 
governed by equation ifT^ , but now (3ru_s — i^Pr (4734 + 
3) - r„/3r„]/[F(4734 + 3) - F„], where 734 = FF„(1 - 
/3r/3r„) is the LF of region 4 relative to region 3. The 
fraction of the decay products deposited into the T-ejecta 
directly, g, is still defined by equations (|16|1 and H17|) . in 
which we take / = min{l, An-f/A}, and f — when 
An-f is negative. 

Energy conservation of the T-ejecta interacting with 
the regions 2 and 4 results in^ 



{U/c^ + M)dT = (7 - 
+ff(F„-F)|dM„| 



r723)(l +lth)dmtr 
- (F„ - F734)dm 

gap 1 



(36) 



° The thermal energy generated in the decay products depositing 
into T-ejecta is small and is therefore ignored here. 
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where U is the thermal energy generated in both the FS 
and RS fronts, which is described by 

dU/c' = (1 - ei) [723(1 + 7th) - l]dmt, 

+ (1 - e2)(734 - l)dmgap, (37) 

dmgup is the differential mass swept by the RS from the 
gap, which can be approximated by 

rfmgap = /Mgap^^^^^, (38) 

where I = 1 for An-r, < A and / = otherwise (see 
equation H17|l for the definition of An-r)- This step 
function is used to characterize whether there is decay 
product falling into the gap region. Mgap is the total 
mass contained in the region (4), which satisfies 

rfMgap = I dMn I ~ dmgap . (39) 

The width of the gap Agap (i.e., the separation of the 
T-ejecta rear and the I-ejecta front), is governed by 

J A - / (^ras - Prj)dR, for An-r < A, , 
rfAgap - j (^j, _ p^^)dR, for An-r > A. ^^^') 

The total mass in the T-ejecta can be calculated 
through 

dM = dnit, + 5|dM„| + dmgap, (41) 

where dm^^ is still defined by equation (|25|) . Notice that 
the expression of dmtr is reduced to dm^^„ when / = is 
satisfied, i.e. there is no decay products deposited onto 
the unshocked medium. This usually does not happen in 
the wind case, but happens sometimes in the ISM case. 

During this stage, the I-ejecta evolves independently. 
Because of the deposition of decay products, the I-ejecta 
is somewhat accelerated. To estimate this effect, we as- 
sume that the I-ejecta moves with a uniform LF Fj. In 
principle, the part that overlaps with the N-ejecta (so 
that deposition of the decay products is possible) should 
move slightly faster than the rest part, and within the 
decay-area different layers may move with slightly dif- 
ferent LFs. But such an effect is unimportant for our 
further discussions about the interaction between the I- 
ejecta and the T-ejecta. Energy and momentum conser- 
vations then yield 

(Mi + t/i/c2)dFi = [(r„-Fi)-Fi(7,ei-l)]/i|dMn|, (42) 

where 

, rn 1 ^gap + Af-r -I- An-f ^ ..„^ 

fi, = max{0,l ^-^ ^ 1, (43) 

and 7i.ci = FiFn(l — PtiPt^) is the relative LF between 
the N-ejecta and the I-ejecta, Ui and Ali are determined 
by dUi/c^ = (7i.ei - l)h\dAQ and dMi = h\dM^\ (Orig- 
inally, Ml = Mj*), respectively. We can then solve Fj at 
any time, which is adopted as the input parameter in the 
next section. 

3.2.2. The collision between the I-ejecta and the 
T-ejecta 

At the beginning of the evolution, one has F > Fj, 
and the I-ejecta lags behind further and further. How- 
ever, the T-ejecta is decelerated by the trail and the ISM 
continuously, so that F decreases with radius. When 



F = Fi the separation between the I-ejecta and the T- 
ejecta is the largest. Later one has F < Fj, and the 
gap between the two ejecta shrinks, until eventually the 
I-ejecta catches up with the T-ejecta. This is a stan- 
dard energy injection process, and a pair of refreshed 
shocks are formed, i.e. a refreshed FS (RFS) expanding 
into the hot T-ejecta, and a refresh RS (RRS) penetrat- 
ing into the I-ejecta. The detailed treatment of such a 
process has been presented before (e.g. Kumar & Pi- 
ran 2000; Zhang & Meszaros 2002). In particular, Zhang 
& Meszaros (2002) considered the interactions between 
three media, i.e. the ISM, the initial fireball shell and an 
injected shell. They considered three emitting shocks, 
including the leading forward shock and the pair of the 
refreshed shocks. Their analysis is pertinent to treat the 
current problem. As shown in Zhang & Meszaros (2002), 
the emission powered by the RFS is unimportant in the 
optical band, (see the short dashed line in their Fig. 4b). 
Since we are mainly focused on the optical lightcurves in 
this paper, that result leads to great simplification of our 
treatment. In our calculations, only the RRS and the FS 
are taken into account. The simplified system includes 5 
parts (see Figure^b) for illustration): (1) the unshocked 
ISM, (2) the unshocked trail (this region may merge to 
region 1), (3) T-ejecta and shocked trail, (4) shocked I- 
ejecta, and (5) unshocked I-ejecta. Notice that we no 
longer separate the unshocked I-ejecta into two regions 
(whether overlap with the N-ejecta), since we treat the 
unshocked I-ejecta as a whole with a same LF Fj. 

Energy conservation of the the system (T-ejecta, the 
shocked trail and the shocked I-ejecta) yields 

([/2/c2 + M2)dF = (7 - F723)(l + 7th)rfmtr 

+ (Fi - F745)(l +7th,i)dmi + .g(Fn - F)|dAf„|, (44) 

where 723 and 745 are the relative LFs between region 
2 and 3, and region 4 and 5, respectively (eq.^Hl), U2 
is the thermal energy generated in the shocks, which is 
determined by 

dU2/c^ = (1 - ei)[(l + 7th)723 - l]dmtr 

+ (l-e2)[(l+7thj)745 -l]dTOi, (45) 
and AI2 is the total mass in the shocked region 

dAh ^ dmtr + dmi + .g|dA/„|. (46) 

Here 7th,i = Ui/A'W, dmi « AfidAi_R/A, dAi_R = 
{Pri—PrTi_s)dR (eq.|221), and Mi is the total I-ejecta mass 
upon collision (according to eq.021)- The velocity of the 
RRS reads /3r„3 « [F/?r(4745 + 3)-Fi/3r:]/[F(4745 + 3)- 

ri]. 

After the reverse shock crosses the I-ejecta, the I-ejecta 
and the T-ejecta merge into one single shell, and the in- 
fiuence of the neutron component essentially disappears. 
The later dynamical evolution of this merged shell is the 
same as that in the neutron-free model, which has been 
studied in many publications (e.g. Huang et al. 2000; 
Panaitescu & Kumar 2002). We do not discuss it further 
in this paper. 

3.2.3. Model Lightcurves 

The dynamical evolution of the system (LF of the T- 
ejecta) and the R-band lightcurve in the ISM-interaction 
case are calculated for x =0i 0.1, 0.5, 1.0, respectively 
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(b) 









Fig. 5. — The early optical afterglow lightcurves of a neutron-fed 
long GRB in the ISM interaction case, (a) The dynamical evolution 
of the region shocked by FS as a function of time. The dotted, dash- 
dotted, dashed, and solid lines represent x = 0.0, 0.1, 0.5, 1.0 
respectively, (b) R-band lightcurves. Line styles are the same as 
in (a). Thick lines represent the total early R-band lightcurves, 
while the thin lines are for RS emission only. For clarity, only the 
RS emission for X = the RRS emission for X = 1 cases are 
plotted. The initial parameters are the same as those listed in the 
caption of Figure|21 except that n = Icm"'' is adopted. 



(see Figure EJ^. In the cases of x = 0-5 and x = 1) 
there is a gap between the I-ejecta and the T-ejecta, and 
a bump is evident in both the F — tobs and F^^ — iobs 
diagrams. This bump is the signature for the I-ejecta and 
T-ejecta colhsion. When x < 0.1, this signature is not 
significant, mainly because the gap is not weU-devefoped. 
The trail is picked up by the I-ejecta continuously, so 
that the treatment is similar to the one to calculate the 
wind case (see H3.1|l . The R-band lightcurve for x = 0.1 
is rather similar to the x = case, except that it is 
relatively brighter at earlier times (e.g. before the RS 
crosses the I-ejecta). The reason is that even in the x = 
0.1 case, the trail is much denser than the ISM, so that 
there are more electrons in the forward shock front that 
contribute to synchrotron radiation. The early R-band 
flux for X = 0.5 and x = 1 is also different from the case 
of X = 0. This is because the early afterglow for a high- 
X case is powered by the T-ejecta, which has a smaller 
mass in the beginning than the I-ejecta. The T-ejecta 
soon enters the deceleration phase, so that F starts to 
drop with time at a very early epoch (see Figure [SJa)). 
In our calculation, it is found out that for high-x case, 
the initial RS that crosses the T-ejecta is too weak for 
any observational consequence. For a comparison, for 
the low-x case (x < 0.1), the more massive I-ejecta is 
decelerated slowly, and it is not significantly decelerated 
within 100 seconds (see Figure|3Ia)). In the meantime, a 
strong RS crosses the I-ejecta, whose emission is included 
in our calculations. This keeps a higher level of the early 
afterglow emission. In all the cases, the rising behavior 
of the early afterglow is because the typical synchrotron 
frequency is above the band. For the high-x case, the 
obvious bump signature is the joint contribution from the 

® Our dynamical evolution of the system for the ISM case is 
also different from that of BOS. The discrepancy is mainly due 
to different assumptions involved. In this work, we assume that 
the I-ejecta moves slower than the N-ejecta, while in BOS, they are 
assumed to be the same. 



FS and the RRS, and latter is the dominant component 
in the optical band. For illustrative purpose, in Figure 
EJb), besides the total lightcurve, we also plotted the 
RRS contribution in the x = 1-0 case, and the initial RS 
contribution in the x = 0.0 case. 

An interesting question is to explore how to use the 
lightcurve to diagnose the existence of the neutron com- 
ponent for the ISM case. For example, the initial rising 
lightcurve of a neutron-free fireball is steep, while for 
a neutron-fed fireball it is much milder because of the 
contribution of the shocked trail emission. The sharp 
bump in the early afterglow stage is a prominent signa- 
ture. However, it suggests that an energy injection event 
happens, not necessarily a proof of the existence of the 
neutron component. Nonetheless, the neutron-fed model 
gives a natural mechanism for post energy injection, and 
it suggests that an early bump is common for long GRBs 
with ISM-interaction as long as x is reasonably large, say 
> 0.5. If future Swift observations reveal that an early 
bump is common, it may suggest that the GRB fireballs 
are baryonic with a large fraction of neutrons. It is worth 
mentioning that for the neutron-free model, an early af- 
terglow bump is also expected due to the transition from 
the RS emission to the FS emission (e.g. Kobayashi & 
Zhang 2003a; Zhang et al. 2003). Nonetheless, the bump 
due to energy injection (presumably powered by the col- 
hsion between the I-ejecta and the T-ejecta) could be 
well-differentiated. For example, the injection bump is 
achromatic, while the FS peak is chromatic, caused by 
crossing of the typical frequency into the band. Finally, 
the RS emission for x = is weaker than the RRS emis- 
sion in the high-x case (FigurelSJ^b)). This is because the 
ejecta is assumed to be cold for the neutron-free case, 
while for the neutron-rich case, it is hot (see equation 
||21)). At the RRS front, the original thermal energy of 
the protons would be also shared by the electrons, lead- 
ing to stronger emission. 

The neutron-fed signatures last only for a few 100 sec- 
onds. Rapid and frequent monitoring (e.g. with ^ 5 s 
exposure interval) of early optical afterglow is needed to 
catch these signatures. Swift UVOT, with an on-target 
time 60-100 s, has the capability to detect such signals. 

4. EARLY OPTICAL AFTERGLOW LIGHTCURVES: SHORT 

GRBS 

Bursts with durations shorter than 2 s may have differ- 
ent physical origin. The leading model for short GRBs 
invokes mergers of two compact objects (e.g., NS-NS 
merger or BH-NS merger, see e.g. Eichler et al. 1989). 
Although there are other suggestions about the origin of 
short GRBs, the merger model was found suitable to in- 
terpret many short GRB properties (e.g. Ruffert et al. 
1997). A pre-concept of such a merger model is that 
the burst site is expected to have a large offset from the 
host galaxy due to asymmetric kicks during the births 
of the NSs, and that the number density of the external 
medium is low, typically O.Olcm^'^. Such a suggestion 
seems to receive support from the afterglow data of the 
recent short, soft GRB 040924 (Fan et al. 2005b). As al- 
ready mentioned in the introduction, in the compact ob- 
ject merger model (especially the NS-NS merger model), 
the outflow is very likely to be neutron-rich (x > 1). 
Here we calculate the early optical afterglow lightcurve 
of a neutron- rich short GRB (x = 1) and compare it with 
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the neutron-free case. 

We take the following typical parameters for a short 
GRB. The burst duration is taken ~ 0.2s, so that the 
width of the N-ejecta is A ~ 6 x lO^cm. The total 
(include both protons and neutrons) isotropic energy is 
i^tot = 2.0 X lO^^ergs. The initial LF of the merged ion 
shell and the neutron shell are Fm = 200, Fn = 300, re- 
spectively The I-eject and the N-ejecta separate at a 
radius « 2F^F2 A/(F2 - F^) « lOi^cm, which is much 
smaller than the typical decay radius ~ 8 x lO^^cm. As a 
result, for short bursts most of the neutron decay prod- 
ucts are deposited onto the ISM (i.e. / ~ 1). This 
is the "clean" scenario discussed by BOS. The result- 
ing C{R) in the trail is very large for R ^ several Rp, 
and a gap still forms between the trail (T-ejecta) and 
the I-ejecta. The calculation is therefore the same as 
the ISM-model for long GRBs (see ^^fov detail). For 
A ^ lO^^cm and Fm ~ 200, the I-ejecta starts to 
spread at i? > lO^^cm (e.g., Piran 1999). This radius 
is much smaller than the deceleration radius, which is 



7x10^"^;, 



16pl/3 ^-l/3p-2/3 



tot, 51' 
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cm for X = 0. In our calcula- 
tion, we take the spreaded width '-^ 5 x lO^^cm when cal- 
culating the RS emission. Also here we perform numer- 
ical calculations as compared with the analytical treat- 
ment presented in Fan et al. (2005b). We found that 
the relative LF between the unshocked I-ejecta and the 
shocked I-ejecta is smaller than the one estimated in Fan 
et al. (2005b), so that the RS emission is further sup- 
pressed. 

The lightcurves are plotted in figure El for x = and 
X = 1. Similar to figure [Sfb), the main difference be- 
tween the two cases is that for the latter the very early 
lightcurve increases more slowly with time. This is due 
to the contribution from the shocked trail. A collision 
also happens as the I-ejecta catches up with the deceler- 
ated T-ejecta. The signature, which is manifested as a 
steeper increasing lightcurve around 20 s, is however not 
as prominent as the long GRB case. There are two rea- 
sons. One is that the T-ejecta and the 1-ejecta now have a 
comparable mass, since essentially all the neutron decay 
products are stored in the T-ejecta (for comparison, for 
long bursts, a good fraction is deposited onto the I-ejecta 
so that the I-ejecta is much more energetic than the T- 
ejecta). Second, the RRS is very weak, and the main 
contribution to the lightcurve is the FS component (for 
comparison, for long bursts, the RRS component domi- 
nates the bump emission). Nonetheless, we identify some 
novel features of a neutron- rich short GRB. The signa- 
ture occurs too early, however, and the global afterglow 
level is very dim due to the small i?tot and n. It is still 
a challenging task to diagnose the neutron component in 
short GRBs. 

5. DISCUSSION AND SUMMARY 

The composition of the outflow that powers GRBs is 
still far from clear. The well-collected late afterglow 
data are unfortunately not suitable to study the GRB 
fireball composition, since the afterglow radiation comes 
from the shocked medium rather than the fireball mate- 
rials. Understanding the fireball composition and hence 

Notice that for a short GRB, as the I-ejecta crosses the shghtly 
decayed neutron shell, the interaction is so weak that there is no 
UV flash predicted for long GRBs (Fan & Wei 2004). 




Fig. 6. — Early R-Band lightcurves for a typical short GRB 
for X = t (solid) and X = (dotted), respectively. The input 
parameters are the same as those listed in the caption of flgure 1^ 
except that n = O.Olcm"^ and Etot = 2 X lO^^ergs are adopted. 



the nature of the explosion requires detailed early after- 
glow data. In general a GRB fireball is composed of 
two distinct components, a baryonic component and a 
Poynting flux component. Current early afterglow data 
suggest that GRB ejecta seem to be magnetized (e.g. 
Fan et al. 2002; Zhang et al. 2003; Kumar & Panaitescu 
2003). Closer modeling suggests that for GRB 990123 
and GRB 021211 in which the reverse shock magnetic 
field is stronger than that in the forward shock region, 
the a parameter, i.e. the Poynting flux to baryonic flux 
ratio, is moderate (Zhang & Kobayashi 2004; Fan et 
al. 2004b). In other two cases of early afterglow (GRB 
021004 and GRB 030418), detailed modehng is needed to 
reveal whether Poynting flux is important. In any case, 
it is likely that the baryonic component is at least non- 
negligible. As discussed in the introduction, the neutron 
component is an inevitable composition within the bary- 
onic component, and the diagnosis of the neutron compo- 
nent is a handle to reveal the significance of the baryonic 
component in GRB fireballs. In the literature, the possi- 
ble neutron signatures discussed include the multi-GeV 
neutrino emission (e.g. Derishev et al. 1999a; Bachall 
& Meszaros 2000; Meszaros & Rees 2000) and the UV 
flash accompanying the 7— ray emission phase (Fan & 
Wei 2004) The detection of these signals are difficult 
due to the limitation of the current detector capability. 

The most prominent neutron signature is likely im- 
printed in the early afterglow. This suggestion has been 
proposed (e.g. Derishev et al. 1999b; B03), but no de- 
tailed calculations have been performed. In this paper, 
we present a first detailed calculation of early optical 
afterglow lightcurves for a neutron-fed fireball. We con- 
sidered both long and short GRBs, and for long GRBs 
we consider both a ISM medium and a wind medium. 
For each model, we calculate the lightcurves for differ- 
ent X (the neutron-to-proton number density ratio) val- 
ues, aiming to study how the neutron component pro- 
gressively change the lightcurve behavior as x increases. 

This prediction may have already been confirmed by the re- 
cent detection of the prompt optical and IR emission (which are 
variable and correlated with the 7— ray emission) during GRB 
041219a (Blake et al. 2005; Vestrand et al. 2005). 
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We confirmed the previous suggestions that the pres- 
ence of neutrons changes the fireball dynamics and the 
lightcurves (Derishev et al. 1999b; B03). Our findings 
are summarized as follows. 

1. For short GRBs, the neutron decay products de- 
posit onto the ISM medium mostly, and the neutron de- 
cay trail clearly separates from the ion ejecta (I-ejecta). 
This is the clean picture delineated in B03. For long 
GRBs, the picture is more complicated. Only part of the 
decay products are deposited onto the medium. The rest 
are deposited onto the Tejecta itself or onto the shocked 
region. The neutron signatures in long GRBs therefore 
require more complicated treatments. 

2. If the medium is a pre-stellar wind (for long 
GRBs), the neutron trail moves slowly (mainly because 
the medium inertia is too large). The trail and the T 
ejecta do not separate from each other, and a forward 
shock propagates into the trail directly. Three compo- 
nents contribute to the final emission, i.e. the forward 
shock, the reverse shock propagates into the I-ejecta, and 
the unshocked trail emission. The latter is significant 
when X is large, since the internal energy of the un- 
shocked trail is large when the medium density is high. A 
typical neutron-rich wind-interaction lightcurve is a char- 
acterized by a prominent early plateau lasting for ~ 100 
s followed by a normal power-law decay (Fig[2J). We also 
show that in the wind case, the IC coohng effect due 
to the overlapping of the initial prompt 7— ray with the 
shocks and the trail (e.g. Beloborodov 2005; Fan et al. 
2005a) suppresses the very early R-band afterglow sig- 
nificantly, and the neutron-fed signature is also dimmed 
(see Figl^b) and Figl^Jc) for a comparison). 

3. If the medium is a constant density ISM (for long 
GRBs) , part of the neutron decay products fall onto the 
medium, and the trail moves faster than the I-ejecta. 
A gap likely forms between the leading trail and the I- 
ejecta. The former forms a distinct trail ejecta (T-ejecta) 
which interacts with the out trail or ISM. The latter 
catches up later and gives rise to a rebrigtening signa- 
ture. Before collision, the radiation is dominated by the 
forward shock emission. During the collision, both the 
forward shock emission and the refreshed shocks (espe- 
cially the refreshed reverse shock) are important. The 
unshocked trail emission is not important in this case. A 
typical neutron-rich ISM-interaction lightcurve is char- 
acterized by a slow initial rising lightcurve followed by 
a prominent bump signature around tens to hundreds of 
seconds (FigEJ. 

4. The picture for short GRBs is similar to the case 
of long GRBs with ISM-interaction. The injection bump 
is less prominent and the refreshed reverse shock is not 
important. A typical neutron-rich short GRB lightcurve 
is characterized by a slow initial rising lightcurve followed 
by a step-like injection signature (FigEJ. 

For all the cases, the predicted signatures can be de- 
tected by the UVOT on board the Swift observatory. 
However, most of these signatures (such as the plateau 
and the bump signature) are not exclusively for neutron 
decay. More detailed modeling and case study are needed 
to verify the existence of the neutron component. 

A neutron-fed fireball involves very complicated pro- 
cesses and some not well-studied physical problems (e.g. 
particle acceleration and emission in the trail) . In princi- 
ple a complicated numerical model is needed to delineate 



the problem. In this paper, we made some approxima- 
tions to simplify the problem and eventually come up 
with a semi-analytical model. Such a treatment nonethe- 
less catches the main novel features of the model. Further 
studies are needed to build up a more realistic model of 
neutron-fed fireballs. 

So far there are four GRBs whose early afterglows are 
detected (Akerlof et al. 1999; Fox et al. 2003a, 2003b; W. 
Li et al. 2003; Rykoff et al. 2004, see Zhang & Meszaros 
2004 for a recent review), i.e. GRB 990123, GRB 021004, 
GRB 021211 and GRB 030418. Modeling GRB 990123 
and GRB 021211 suggests that these two bursts are born 
in an ISM medium (e.g. Panaitescu & Kumar 2002; Ku- 
mar & Panaitescu 2003). The early afterglow of these 
two bursts are dominated by fast-decaying optical flash, 
usually interpreted as the reverse shock emission compo- 
nent. No neutron-signature discussed in this paper has 
been discovered. There are two possible reasons for this. 
First, the ejecta of these two bursts are likely magnetized 
(Fan et al. 2002, 2004b; Zhang et al. 2003; Zhang & 
Kobayashi 2004; Kumar & Panaitescu 2003; Panaitescu 
& Kumar 2004; McMahon et al. 2004). This would dilute 
the neutron signature discussed in this paper (in which 
(7 = has been assumed throughout). Second, modeling 
suggests that the initial LF of the I-ejecta of these two 
bursts are large, around 1000 (e.g. Wang et al. 2000; 
Soderberg & Ramirez-Ruiz 2002; Wei 2003; Kumar & 
Panaitescu 2003). Since the N-ejecta could not be ac- 
celerated to a very high LF (Bahcall & Meszaros 2000), 
the N-ejccta likely lags behind the I-ejecta, and its in- 
fluence on the fireball dynamics is minimized. The early 
afterglows of other two bursts GRB 021004 and GRB 
030418 are not easy to interpret within the standard re- 
verse shock picture. One suggestion is that both bursts 
are born in a wind environment (e.g. Li & Chevalier 
2003; Rykoff et al. 2004). If this is the case, the available 
earliest detection at ^ 300s is too late to catch the trail 
emission predicted in this paper. Alternatively, they are 
potentially interpreted in an ISM model with the neutron 
signature discussed in this paper. With the combination 
of the injection bump and the forward shock bump (e.g. 
Kobayashi & Zhang 2003a), one might be able to achieve 
a broad early bump to interpret the early afterglow be- 
havior of GRB 021004 and GRB 030418. Finally, these 
two bursts may be modeled by a high-cr flow (e.g. Zhang 
& Kobayashi 2004 for more discussions). More detailed 
modeling is needed to verify these suggestions. 

With the presence of magnetic fields, the acceleration 
and interaction of the neutron component may have some 
novel features. In the magnetization scenario, a two com- 
ponent jet is likely, with the wide less coUimated jet be- 
ing due to the mildly relativistic neutrons (Vlahakis et 
al. 2003; Peng et al. 2004). Peng et al. (2004) suggest 
that the decay products in the wider neutron compo- 
nent would give rise to an injection bump in the after- 
glow for an observer on-beam the narrow core jet. In 
such a picture, for an observer off-beam of the narrow 
jet but on-beam the wide jet, the observer would see an 
orphan-afterglow- type event. Since the LF of the wide 
component is about 15, the neutrons would decay at a 
typical distance of ~ 4 x 10^'' cm. The decay products 
shock into the ambient medium and get decelerated at a 
further distance (the deceleration distance). For some- 
what optimistic parameters (cc = 0.3, eb = 0.1, p = 2.3, 
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n — I and = 1/3), the forward shock typical syn- 
chrotron frequency at the typical neutron decay radius 
is Vm ~ eV for an ISM medium, or i/m ~ 0.1 keV for 
a wind-medium. In the ISM case, the optical lightcurve 
increases rapidly with time and reaches a peak at the de- 
celeration time (typically hours after the burst trigger). 
In the wind case, since the soft X-ray band is above both 
Vrn and Vc, the lightcurve between the typical decay time 
and the deceleration time is essentially flat (oc t(^~^'^/^), 
giving rise to a flat soft X-ray plateau lasting for hours. 
This would be an interesting signature for the neutron- 
rich two-component jet model in the wind environment. 
However, without a 7-ray or hard X-ray precursor, de- 
tecting such a signature is challenging. 
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APPENDIX 

APPENDIX: CALCULATIONS OF SYNCHROTRON RADIATION 

Here we present the standard treatment of the synchrotron radiation of electrons with a power-law energy distribution 
(e.g. Sari, Piran & Narayan 1998, and Piran 1999 for a review). This is valid for all the cases invoking shocks, and 
it is also valid for one extreme model of trail emission (as discussed in Appendix^. As usual, we introduce the 
equipartition parameters Ce — 0.1 and es — 0.01 as the energy fraction of the electrons and magnetic fields in the 
local thermal energy in the energy dissipation region (e.g. shock front), respectively. The electrons are assumed to be 
distributed as a power law, i.e. dNc/djc « 7(7^ for 7c > 7c, m, where 

7e,m = ec(mp/mc)[(p - 2)/{p - 1)] [7^(1 + 7th) - 1] + 1, (Al) 

TUp and TTic are the rest mass of proton and electron respectively, p ~ 2.3 is the typical power law index of the electrons, 
7r is the relative LF between the shock upstream and downstream, and 7th is the average random LF of protons in the 
upstream. When deriving this expression, we have taken the downstream random LF as 7th, d = 7r(l + 7th,u) — 1 (where 
7th, u = 7th in the above expression), which is approximately valid when 7th. d 3> 7th, u- The exact relativistic shock jump 
conditions for a hot upstream (Kumar & Piran 2000; Zhang & Meszaros 2002) leads to 7th, d + 1 — (4/3)7r(l -I- 7th, u) 
when 7th, d ^ 7th, u is satisfied. So our treatment is valid to order of magnitude. 
The co-moving magnetic field strength in the shock front can be estimated by 



B 



' ~ y^327reB7r[7r(l + 7th) - l]n^m.pc'^ 



(A2) 



where riu is the co-moving proton number density in upstream. At a particular time, there is a critical LF of the 
electrons, 

BtTTOcC 



7e 



{l + Y)aTTB'^t 



(A3) 



above which the electrons are cooled, F is the bulk LF of the shocked region. Here ctt is the Thomson cross section, 
Y is the Compton parameter (e.g., Wei & Lu 1998, 2000; Panaitescu & Kumar 2000; Sari & Esin 2001; Zhang & 



Meszaros 2001) which reads Y 



-, where x is the radiation coefficient of the electrons, so that the total 



radiation efficiency is eiad = xe^, where x — min{l, (1 + fc)(7c.m/7c,c)'''^^''/[^(3 ^p)]}i where A: = 3, 1 are for the ISM 
and wind respectively. This last equation could be derived as follows. 

Assuming that the outflow is in the slow cooling phase, one has dNc/d^jc ~ Ao(7c/7e,m) for 7c, m < 7c < 7c. c 
and dN^/d'jc = A^o7c,c7c~m(7c/7c,m)"'P+^' for 7c,c < 7c < 7c,m, where 7c,m 10^(B7lG)-i/2 ^ ^^^^^ xhus the total 
emitting power (including the IC component) satisfies 



Ptot = (4/3)(l + r)7VorVT(S'V8^)c7cP„J / 7c~'''"'^t^7c + 7c 

« (B'V6^)(1 + Y)N,T^a^nl^^^l-^/[{p - 2) (3 - p)]. 
On the other hand, the total energy of the electrons reads 

/■7c, c /•7c, M 

W = 7Vor7P,.„™cc2 [ / 7c"'''~'^rf7c + 7c,c / l-^d-f, 

•'7o,m •'7o,c 

^NoT-il^m,^/{p-2). 
The ratio of the fresh electrons {dNc) to the total electrons (A'c) satisfies 

dN^/N^^ kdt/[{l + k)t], 



-(p-i) 



d-fc 



(A4) 



(A5) 
(A6) 
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where fc = 3, 1 are for the ISM and wind respectively. In deriving equation ljA6p . relations oc R^, dR ~ AT^cdt 
and R ~ 4(1 + k)T^ct have been used, and the outflow is assumed to be in the adiabatic phase. So the total energy of 
the fresh electrons is dW « WdNc/Nc, the corresponding total input power (including that of protons) is 



Pfrcsh = dW/{e,dt) » fciVor7,',^mec7[ee(l + k){p - 2)t], 
Combining with equation ljA3p . wc have 

X = -^tot ^ (1 +fc) 7can .p-2 
~ eoPfrcsh fc(3 - p) 7c,c 

This more precise result is slightly larger than that presented in Sari & Esin (2001). 
The typical synchrotron radiation frequency and the cooling frequency are estimated by 



(A7) 



(A8) 



and 



2(1 + z)TrmeC 



2(1 + z)7rmoc' 



(A9) 



(AID) 



respectively, where e is the electron charge, z is the redshift of the GRB. 

Another potentially important frequency involved (in the wind case) is the self- absorption frequency ly^. Here we 
follow Rybicki & Lightman (1979, p. 188-190) for a standard treatment: In the co-moving frame of the emission region 
(all physical parameters measured in co-moving frame are denoted with a prime) , the absorption coefficient a'^^, scales 
as a'^, (X ty'~(P+^^/^ for i/^ = max(t'^j, i/'^) > i/' > v'^ = min(^'^, i/^), and as a'^, cx v'~^^^ for v' < ly'^. For i/^^ > > v^, 
one has 



while for ly' < v'^, one has 



Here A 



( 3e \ 



p/2 



{m.c^ f-^KB'^^+^y^T (^) r (^), where K : 



2(p-l)n^„7P-V3, n. 



(All) 

(A12) 
is the down- 



SttttIc \2-!rm'ic^ J V-c^ I \ 12 ) " \ 12 7' -v/' ^J"-y^lp / '^i '"70 

stream proper proton number density. 7p = min(7o,m, 7c, c), r(x) is the Gamma function. The above general treatment 
is valid when we take p = p for the slow cooling and p = 2 for the fast cooling. The self-absorption optical depth can 
be calculated by 



-(^') 



a',dR'. 



(A13) 



When considering shocks, we can approximate / a'^,dR' w a'^,AR' , where AR' is the width of the shocked region, so 
that n^^AR' ~ Nc/4:TtR'^, where Nc is the total number of emitting electrons. The co-moving self-absorption frequency 
v'^ can be derived from r(z/^^) = 1, and the observed absorption frequency is z^a = Tv'^/{1 + z). 

The synchrotron flux as a function of observer frequency can be approximated as a four-segment broken power law. 
For < minji/c, J^m}, see Eqs. (4-5) of Zhang & Meszaros (2004). For > minji^c, J^m}, one has 



F « F 



for the fast cooling case, and 



(^'c/^^a)^(^^/^'c)^ (i^ < t^c), 

(^/a/^^c)-l/2(^^/^^a)^/^ (i^c < ly < Z^a), 

{vJv.r^^\ylyu.)-^^\ {y > y.n). 



(A14) 



F ^ F 



(^)-^(^)5/2^ 



(y > i/c), 



(A15) 



for the slow cooling case. Here F^r 



3\/3^p(l+z)jV 



TB, <l>p is a function of p (e.g. for p ~ 2.3, $p ~ 0.6) 



(Wijers & Galama 1999), and Z3l being the luminosity distance. 
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APPENDIX: TRAIL EMISSION 

One novel feature of the neutron-fed GRB afterglow, in particular in the wind-medium case, is the emission from the 
neutron decay trail itself (without being shocked) . This is also one of the main challenges faced in preparing this work. 
For example, it is unclear (1) what is the mechanism of trail emission, (2) how much thermal energy is distributed to 
electrons and magnetic fields, and (3) how the emitting electrons are distributed. 

The first two uncertainties also apply to the shock acceleration case. Our treatment here closely follows the stan- 
dard shock model. We first assume that the trail emission mechanism is also synchrotron radiation, and assign two 
equipartition parameters, i.e. ~ 0.1 and ee 0.01 for the electrons and the magnetic fields, respectively. For the 
third uncertainty, while there is a standard paradigm (i.e. Fermi acceleration) in the shock case so that a power-law 
distribution of electron energy is justified, it is unclear whether it is the case in the trail. Rather than speculating the 
possible electron distributions, in this paper we consider two extreme situations, which may bracket the more realistic 
electron distribution case. In the first model, we still consider a pure power-law model that is completely analogous 
to the shock case, i.e. for an average electron LF we get the minimum electron LF ^e,ra = [(p — 2)/(p — 1)]7g, 
and dNe/dje cx (7e/7e,m)~^ for 7e > 7e,m- In the second model, the electrons are assumed to be mono-energetic, 
i.e. all electrons have a LF oi %. A more realistic version of the mono-energetic model is the relativistic Maxwellian 
distribution model, i.e. dN^/d'^c oc 7^ exp(— 7cTOoC^/fc^c)j where Tc is the temperature of the plasma. For an observer 
frequency far above several xj/e, where Vc = l{R)lc^B' /2ti{\ + z)mcC, the observed flux for both the mono-energetic 
model and the Maxwellian model is much dimmer than the one in the power law model, where "f{R) is the bulk LF of 
the trail. Below i^e, on the other hand, the synchrotron emission properties are quite similar for various distribution 
models. Fortunately, for the wind-interaction case when the trail emission becomes important, Pc is indeed above the 
R-band frequency z/r — 4.6 x lO^^'lIz, so that the uncertainty to calculate the trail emission lightcurve is small (see 
Fig|3. 

The magnetic field generated in the neutron-front can be estimated by 

-Bjr ~ ^87reB7th?itrmpc2, (Bl) 

where ntr has been shown in equation H12|l . Please notice the different expression for B' and We also assume 
ee ~ 0.1 in the trail, which may be conservative. Even with this estimate, the trail emission is found to be strong 
enough to be detectable in the wind case (see Figl^for detail). 

A self-consistent calculation for the trail emission is quite complicated. In this work we make the following approx- 
imations, (i) After the medium is "ignited" by the neutron front at 

tig„ - R/2Tlc, (B2) 

the trail located at R (but moving with 7(i?)) continually contributes to the observed flux until it is "terminated" by 
the I-ejecta shock front, (ii) The trail is divided into many sub-layers, each moving with 7(i?) without interaction. 
The observed trail emission is a sum of the independent radiation from these sub-layers. During the "lifetime" of each 
sub-layer, the thermal LF of a particular electron cools as 7c(i) = 7c,o/[l + (1 + Y){aTB[j.'^ /QTTmcc)jc,oj{R){t — ^ign)], 
where 7e,o is the initial LF of the electron. For the pow er la w mo del, the electrons are assumed to be distributed as a 
broken power law (considering cooling), but equations lAip and (|A3p are now replaced by 

7c,m(t) - 7c,m/[l + (1 + Y){<jtB'J /67rm,c)-f,^n,j{R){t - Ugn)] (B3) 

and 

7e,c(t) = 67rmec/(l + Y)[aTB'J-f{R){t - Ug^)], (B4) 

respectively. 

In the ISM case, for u-^sm — 1 cm~^ (the favored value for the current multi- wavelength afterglow modeling), the 
trail emission never becomes dominant in the early R-band afterglows. This is contrary to the wind case, in which 
C{R) is only moderate while 7th S> 1, even at i? <C i?/3. The synchrotron radiation from the trail is therefore strong, 
which contributes significantly to the early R-band afterglow lightcurves (see FiglS)). 
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